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Effects of Boundary Conditions and Initial 
Out-of-Roundness on the Strength of 
Thin-Walled Cylinders Subject to 


External Hydrostatic Pressure 


By G. D. GALLETLY! anv R. BART? 


Using classical small-deflection theory, an investigation 
was made of the effects of boundary conditions and initial 
out-of-roundness on the strength of cylinders subject to 
external-hydrostatic pressure. The equations developed 
in this paper for initially out-of-round cylinders with 
clamped ends, and a slightly modified form of the equa- 
tions previously derived by Bodner and Berks for simply 
supported ends, were applied to some actual test results ob- 
tained from nine steel cylinders which had been subjected 
to external hydrostatic pressure. Three semiempirical 
methods for determining the initial out-of-roundness of 
the cylinders also were investigated and these are described 
in the paper. The investigation indicates that if the 
initial out-of-roundness is determined in a manner similar 
to that suggested by Holt then the correlation between 
the experimental and theoretical results is quite good. 
The investigation also indicates that while the difference 
in collapse pressures for clamped-end and simply sup- 
ported perfect cylinders may be quite considerable, this 
does not appear to be the case when initial out-of-round- 
nesses of a practical magnitude are considered. 


INTRODUCTION 


EVERAL analyses have appeared in the literature for the 
S elastic buckling of a thin cylindrical shell subject to external 

hydrostatic pressure (1). The majority of these analyses 
have been based upon the classical small-deformation theory of 
thin shells and have assumed a geometrically perfect, stress-free 
structure prior to loading. ‘The correlation obtained between 
these theories and experimental results has been good for long 
cylinders but rather poor for short cylinders. Efforts are cur- 
rently being made by several investigators to explain the dis- 
crepancy between theory and experiment in the short-cylinder 
range by the use of large-deflection theory. At this date, how- 
ever, it is not known by how much this discrepancy will be re- 
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duced as all the final reports on their work have not been pub- 
lished. 

One possible cause for the discrepancy between theoretical and 
experimental results can be ascribed to the initial out-of-round- 
ness of the cylinders, and a number of investigations, using small- 
deflection theory, already have been made on the effect of initial 
irregularities on the collapse pressure of cylinders subject to ex- 
ternal hydrostatic pressure (2, 3, 4). As was to be expected, these 
analyses showed that the initial irregularities reduced the failure 
pressures below those of the perfect cylinders. However, when 
these analyses were applied to some models which had been tested 
experimentally, they predicted failure pressures which were less 
than three quarters of those observed experimentally. Since 
these analyses had assumed simple supports at the ends of the 
cylinders and it was probable that the boundary conditions of the 
models were somewhere between the extremes of simple supports 
and clamped ends, it was of interest to investigate the reduction 
in collapse pressure of clamped-end cylinders due to initial irregu- 
larities, to see if the assumed boundary conditions would explain 
the discrepancy between experiment and theory. Also, the 
analyses assume that the initial out-of-roundness in the cylinders 
is similar to one of the modes into which a perfectly circular 
cylinder of the same dimensions would buckle, and actual shells 
never satisfy this condition. It thus seemed desirable to in- 
vestigate the various simplified methods that have been sug- 
gested for determining the initial out-of-roundness of the cylinders 
to see what effect these had upon the computed failure pressure. 
These methods are described in the paper. 

It is also of interest to note that there are other limitations in 
the existing linearized theories. These are: 


(a) The fact that in these problems there usually exist other 
buckling pressures close to the minimum buckling pressure. 
Thus, use of only that mode of initial out-of-roundness which 
corresponds to the minimum buckling pressure is really only de- 
fensible at pressures very close to the minimum buckling pressure. 

(b) The simple yield criterion used to predict failure. 


This point is discussed under the section entitied ‘‘Assumptions 
Made in Analysis.” 

The authors have not investigated the foregoing factors but 
hope to do so in the future. 

The approach used in this paper is similar to that of Bodner and 
Berks (3), except that instead of using a Donnell-type equation 
Galerkin’s method was employed in conjunction with a modified 
Donnell-type equation. The initial out-of-roundness pattern 
assumed by Bodner and Berks was of the form 


; 4x 
wo = e sin m 8 cos — 


(origin at mid-length) while that assumed by the present authors 
was 
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(origin at one end of the cylinder). Thus, in both cases, the 
initial out-of-roundness satisfied the relevant boundary condi- 
tions and was also similar in form to one of the assumed buckling 
modes. The magnitude of the initial out-of-roundness at the ends 
and center of the cylinder was also the same in both cases, 

The Bodner-Berks solution and that presented herein thus 
represent lower and upper bounds for the effect of initial eccen- 
tricities on the collapse pressures of elastically supported un- 
stiffened cylinders, when the initial eccentricities have the same 
shape as one of the assumed buckling modes of the perfect 
cylinder. While the two solutions are not exact, they should 
provide good approximations to the exact solutions, One of the 
limitations of using Donnell’s equation is that the number of cir- 
cumferential lobes should be fairly high, and thus the results will 
be slightly in error for very long cylinders which buckle into two 
or three circumferential lobes. 

The final results of the investigation are given in Fig. 2 and in 
Tables 2, 4, and 5. It can be seen that the correlation between 


experiment and theory is quite good when method (c) is used to, 


determine the initial out-of-roundness of the models. [Method 
(c) is similar to that suggested by Holt (5).] However, it is not 
claimed that the results give a complete answer to the problem 
and more work of both an experimental and theoretical nature is 
required. 


METHOD oF ANALYSIS 


The modifications to Donnell’s equation brought about by 
initial eccentricities in the shell have been presented by Bodner 
and Berks (3), and, prior to them, by Cicala (6). The equations 
also have been derived by the authors in the Appendix, using a 
somewhat different approach to that adopted by the previously 
mentioned authors. For the case of uniform external hydrostatic 
pressure applied on all sides of an imperfect cylinder the relevant 
equations are, from Equations [19a], [19d], [21], and [24] in the 
Appendix 
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= thickness of shell 

= mean radius of shell 

= applied hydrostatic pressure 

= modulus of elasticity 

= Eh*/12(1 — y?) 

= Poisson’s ratio 

= stress function of the total membrane stresses 
°2 1 °? 2 
=— SS, . 4 a es 

(= eae UES eRe 

initial radial out-of-roundness (+ inwards) 


u,v, = elastic axial, tangential, and radial (+ inwards) 
displacements of the imperfect cylinder minus the 
uniform compression experienced by a perfect 
cylinder (see Equation [18] in the Appendix), 


The subscripts x and @ indicate partial differentiations with re- 
spect to those variables. 

The patterns assumed for w and wo, and which satisfy the 
boundary conditions for clamped-end cylinders, were as follows 


2 
w= B sin md | 1 — cos | 


half amplitude of w-displacement 

maximum value of initial radial out-of-roundness 
number of circumferential waves 

unsupported length of shell 

axial and angular co-ordinates 
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If the Expressions [2] happen to be an exact solution of the prob- 
lem, then they will satisfy the differential equation of equilibrium, 
Equation [la], exactly. However, as both w and Wo were chosen 
to satisfy the boundary conditions rather than the equilibrium 
equation, this, in general, will not be the case. The resulting ex- 
pression will be a function of z and @ which we shall denote by Q. 
Galerkin’s equation for determining the relations between the co- 
efficients B and e is then 


2x Ll 
if f Q sin 20 [: — cos “| Rdédz = 0..... [3] 
ty) 0 L 


where 7 assumes the values 1,2,3...... 

For 1 ~ m Equation [3] will be found to be zero identically. 
For 7 = m the following relation between B and e, obtained from 
Equation [3], will be found to hold 


where por is given by the expression 
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The smallest value of the buckling pressure of the perfect 
cylinder por is found by minimizing Equation [5] with respect to 
m. A relation similar to that expressed by Equation [5] has re- 
cently been presented by Nash (7) using an energy method. A re- 
lation similar to: Equation [4] was also obtained by Bodner and 
Berks for simply supported imperfect cylinders. 

Thus, from Equations [2] and [4], we obtain the following ex- 
pression for w 


The bending moments in the shell can then be calculated from 
the relations 
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The maximum bending stresses are then given by 


%% ™ 


6 6 
+ a (Ms)max} Cg (Mo)max 
To obtain the total normal stresses we now have to add the 
membrane stresses to Equation [8]. To determine these latter we 
solve Equations [1d] and [6] for the stress function F (periodic 


terms only). The total membrane stresses are then given by 
pR Fog pk 
a a =—— + F,,.... (9 
Ooms 2h Ri’ Tmd h a Wy [ ] 


The total normal stresses are obtained by adding algebraically 
Equations [8] and [9]. The greatest normal stresses occur at 
mid-length of the cylinder (cr = L/2) and where sin m0 = +1 
(trough and crest points of the lobes), At these points the twist- 
ing moment M,, is zero and thus the normal stresses are principal 
stresses. The absolute maximum normal stresses occur at the 
outer shell wall for the trough points. 

Having obtained the maximum principal stresses o, and g¢ in 
terms of the initial out-of-roundness, the geometric parameters 
and the applied external hydrostatic pressure, we now employ the 
octahedral shear-stress criterion of failure (which gives the same 
results as the Hencky-von Mises criterion of failure), viz 

Oe 3 OP oe GP Ulsan gccconoc000 
where a, is the yield point of the material. Substitution of the 
maximum principal stresses 7, and og in Equation [10] then gives 
an equation relating the initial out-of-roundness, the geometric 
parameters of the shell, the yield point of the matexzial and the 
pressure at which the shell begins to yield p,. 


It should be noted that instead of using the yielding criterion 
given by Equation [10] where o, and g¢ are principal stresses, it 
is more accurate to use the expression 


oy? = o9? + 0,2 — dec, + 37,93 


where now @,, 9, and 7,9 are the normal and shear stresses at any 
point and which are functions of zand 6. Yielding will first occur 
in the cylinder for those values of z and 6 which maximize the 
right-hand side of Equation [10a]. However, to compute these 
values of z and @ by differentiation of equation [10a] involves 
more complications than seem warranted. Trials indicate that the 
stress condition at the outer sbell wall for trough points of a lobe is 
probably as unfavorable as anywhere else. As mentioned earlier 
the twisting moment M,¢ is zero at these points and thus Equa- 
tion [10a] reduces to Equation [10]. 

As we shall later present curves of p,, the pressure at which shell 
yielding commences, versus ¢/h, the initial eccentricity- 
shell thickness ratio, for both simply supported and clamped-end 
cylinders, we have summarized the results obtained in this paper 
and those obtained by Bodner and Berks in Tablel. (We have 
added a few terms to the latter solution, as Bodner and Berks 
neglected the periodic terms in Equation [9]). Also, we used 
the expression we/R? for the circumferential change in curva- 
ture instead of (weg/R*? + ve/R) which was used by Bodner 
and Berks when computing the bending stresses due to initial 
out-of-roundness. Our expression is consistent with the expres- 
sion used in deriving the approximate equations of equilibrium 
for an initially out-of-round cylinder (see Appendix) and also is the 
same as that used by Donnell in his work on perfect cylinders (8). 
The effect of using wgg/R? instead of (we9/R? + vg/R) is to 
eliminate the quantity f which appears in the equations de- 
veloped by Bodner and Berks.) It can be seen that the final 
equation (Equation [12] in Table 1) relating the initial out-of- 
roundness and the pressure to cause first yielding is essentially 
the same for both clamped-end and simply supported cylinders, 
except for the factor of 2 required by the definition of e as the 
maximum initial out-of-roundness and the slight changes in 
definition of the quantities K and 8 appearing in that equation. 
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Assumptions Mapw IN ANALYsIS 


Besides the approximate strain and change in curvature-dis- 
placement relations used in the derivation of the Donnell equa- 
tion for an initially out-of-round cylinder, .and the assumption 
that it is initially stress free, it might be useful to list some of the 
other assumptions that have been made both in the analysis of 
this paper and that of Bodner and Berks. These are as follows: 


1 That the circumferential membrane stress in the shell is 
constant along its length and equal to —pR/h, whereas this is 
actually not the case. 

2 The assumption that the initial out-of-roundness is small 
(with an order of magnitude of one shell thickness) is symmetric 
with respect to the center line of the shell, and has the same form 
as one of the buckling modes of a perfect cylinder with the same 
shell dimensions. Actual shells rarely, if ever, satisfy the last two 
requirements and so the question arises as to how the initial out- 
of-roundness should be measured. 

3 The assumption that failure occurs (appearance of visible 
lobes) when the most highly stressed points in the cylinder start 
to-yield. Actually, failure does not occur until plastic regions 
form at the trough and crest points of the lobes. The pressure re- 
quired to produce these yield zones is greater than that at which 
the most highly stressed points begin to yield and neglect of this 
effect therefore underestimates the strength of the shells. An 
adequate theory to take this effect into account has not been de- 
veloped as yet, but, as for beams, presumably the ratio of the 
pressure to cause first yielding to the pressure required for the 
formation of plastic regions depends on the relative magnitudes 
of the direct stresses in the cylinder wail and the bending stresses 
resulting from initial out-of-roundness. 

4 The assumption that Poisson’s ratio is a constant and equals 
0.3. 

As it is intended to apply the analyses developed for unstiffened 
cylinders to stiffened cylinders which failed by buckling between 
ring stiffeners, it would be well if we enumerated the additional 
assumptions that were made. These are: 

5 That the stiffening rings at the ends of any bay are per- 
fectly circular; i.e., they do not have any initial out-of-roundness. 
This never occurs in practice, of course, but should not be too 
serious if the ciroularity of the stiffening rings is very much better 
than that of the shell, or if the predominant mode of initial out-of- 
roundness in the rings is very different from the predominant 
mode in the shell. 

6 As for unstiffened cylinders, the circumferential membrane 
stress in the perfect cylinder is assumed to be —pR/h, whereas 
it actually varies along the length of the shell. A more correct 
representation of the stress distribution would-be obtained by 
using the analysis of von Sanden and Ginther (9), or more ac- 
curately etill, that of Salerno and Pulos (10). 


METHops oF DETERMINING INITIAL OuT-or-ROUNDNESS 


As mentioned hitherto, the analyses assume that the initial 
out-of-roundness is symmetrical about the mid-length of the shell 
and that ite circumferential variation is in the shape of one of the 
buckling modes of a perfect cylinder. Actual shells do not meet 
either of these requirements. If we do not make a harmonic 
analysis of the initial out-of-roundness, and also extend the theory 
to account for the various harmonic components, the question 
arises as to how we shall measure the quantity e, which is defined 
in the analysis. as the maximum initial out-of-roundness when its 
shape is similar to one of the buckling modes. As far as the 
authors are aware three simplified, semiempirical methods for de- 
termining the initial out-of-roundness have been proposed in the 
literature so far. These are (see Fig. 1): 

(a) The centroid of the initial circularity contour ia first deter- 
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mined. Then the angle x/m, where m is the number of lobes into 
which the perfect cylinder would buckle, is calculated. A sector 
of a circle, subtending this angle a /m is then drawn on transparent 
paper and placed with its apex at the centroid of the initial circu- 
larity contour. The sector is then rotated so that it traverses the 
entire circumference of the circularity contour until the location is 
found at which the maximum difference between the two sector 
radii occurs. The initial out-of-roundness is then taken as this 
maximum difference. This method of determining the initial 
out-of-roundness is essentially that proposed by Saunders, Tril- 
ling, and Windenburg (11, 12). 

(b) Both the centroid and the area of the initial circularity con- 
tour are determined. The radius R, of the circle whose 
area is the same as that of the initial circularity contour 
is then determined. A circle, with center at the centroid of 
the initial circularity contour and of radius R,,, is then drawn. 
The initial out-of-roundneas is then taken as the maximum value 
of [R, — R,,] and [R,, — R;], where R, and R, are the radius vec- 
tors from the centroid to points on the initial circularity contour 
which are exterior and interior, respectively, to the circle of radius 
R,,. A method for determining the initial out-of-roundness simi- 
lar to the foregoing has been suggested, among others, by Bodner 
and Berks (3). 

(c) As in (b) both the centroid of the initial circularity contour 
and the radiys R,, of the mean circle are determined. Also, as in 
(a), the angle 7/m is calculated. The arc length of one half-lobe 
is then obtained as(7/m)R,,. This arc is then moved around the 
initial circularity contour with its end points always in contact 
with the contour. The initial out-of-roundness is then taken as 
the maximum radial distance between the circularity contour and 
thearc, This method for determining the initial out-of-roundness 
is somewhat similar to the method proposed by Holt (5). 


NuMERIcaL REesutts 


In this section we present the results obtained by applying the 
analysis developed for simply supported imperfect cylinders (3), 
and its extension to clamped ends, to nine steel welded cylinders 
that have been tested at the Taylor Model Basin. At failure, all 
these models had lobes which partially covered the circumference 
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TABLE 2 GEOMETRIC RATIOS, YIELD POINTS, EXPERIMENTAL AND THEORETICAL COLLAPSE PRESSURES FOR STEEL CYLINDERS 


(1) (2) (8) (4) (5) (6) 
Expt. 
buokling 
Modal us cy Number of pressure, 
no. 2R 2R (psi) stiffeners psi 
50 0.125 0.0028 39000 2 internal 175 
83 0.250 0.0032 31000 2 internal 133 
66 0.500 0.0048 40000 none 2356 
71 0.500 0.0065 44000 none 3276 
42 1.000 0.0033 48000 none 58 
61 2.000 0.0040 39000 none 48d 
BR-1 0.184 0.0024 61700 5 external 
BR-4 0.250 0.0049 50600 5 external 390 
BR-5 0.184 0.0023 54400 6 external 80 


«© Numbers in parentheses are the number of lobes at failure. 


(1) (8) (9) (10) 
Pressure to 
Maximum Theoretical clastic buckling cause axisym- 
sustaining pressure from Eq. [11] of metric yielding 
pressure, ——Tablell\\psi: of stiffened 
pei Simple supports Clamped ends oylinders ° 
195 (18, 19)¢ 301 (17) 451.5 (18) 264 
189 (138, 14) 192.5 12) 262 50) 225 
235 (8, 9) 283 8) 395 ae 6a0 
827 (?) 586 8) 804 9) 
59 (8) 51.2 7) 73.4 8) 
48 (65, 6) 41.1 5) 60.2 6) ae 
107 1) 135 15) 186 res} 322 
390 (10) 613 11) 848 12 560 
95 (14) 120 15) 166 17) 284 


» No distinction made in Windenburg’s notes between buckling and maximum sustaining pressures for these models; henve, they were assumed to be equal 


in the above table. i i‘ 
© Caloulated for median surface of shell at mid-bay length using analysis of 


TABLE 8 MAXIMUM ¢/h-VALUES OBTAINED BY DIFFERENT 


von Sanden and Ginther in conjunotion with octahedral shear stress oriterion 


METHODS FOR DETERMINING INITIAL OUT-OF-ROUNDNESS 


1 (2) (8) (4) (5) (8) 
(1) $$$ $—$$____—__—_———Method (0) —_——_—_,, 
Method (a)—— ——(i) Outward———. __ ———(ii) Inward——— -—Marx. of (i) and a 
Model no C.E. 8.8. Method (b) C.E. 8.8. C.E, 8.8. C.E. 8.8. 
50 0.186¢ 0.140 0.175 0.016 0.017 0.074 0.078 0.074 0.078 
33 0.097 0.104 0.008 0.067 .078 0.085 0.093 0.085 0.003 
66 0.043 0.047 0.041 0.041 0,045 0.032 0.037 0.041 0.045 
71 0.1389 0,143 0.118 0.008 0.011 0.090 0.008 0;090 0.098 
42 0.267 0.272 0.191 0.080 0.091 0.114 0.121 0.114 0.121 
61 0.125 0.180 0.082 0.027 0.083 0.022 0.027 0.027 .083 
BR-1 (6, Sis 0.575(6)> 0.615(8) 0. 554(2) Parana} 0.375(6) 0. 267(6) Guegerg) 0.819 8 0.375 R 
BR-4 [4, 10] 0.164(4) 0.176(4) 0.128(2) 0.063(4) Gece 0.073(4) 0.082(4 0.073(4 0.082(4 
BR-5 [4, 6, 10] 0.741(4) 0.758(4) 0.502(4) 0. 236(10) 0. 285(10) 0.210(4) 0, 247(4) 0.236(10) 0. 285(10) 


® Tabulated values accurate to approximately +5 per cent. 
b Numbers in square brackets in 


cate the stations at which lobes first sppeared in the multibay stiffened oylinders. Numbers in parentheses are the sta- 


tions at which the maximum ¢/h-values occurred, according to the method used. 


TABLE 4 a OF PREDICTED PRESSURES 


P8I) FOR OCCURRENCE OF A VISIBLE LOBE—USING THE DIFFERENT 


METHODS FOR DETERMINING INITIAL OUT-OF-ROUNDNESS—WITH THE EXPERIMENTAL PRESSURES (PSI) 


M) (2) at (4) (5) (8) (7) 8) 
Model —Method (a)— —Method (b)—~ Method (c)—- —Col. (2)/Col. (5)— —Col. (8)/Col. (5)— —Col. (4)/Col. (6)— 
no. C.E. 8. C.E. 8.8. C.E. 8.8. Expt. C.E. 8.8. C.E. 8.8. C.E. 8.8. 
50 351 134 137 125 180 160 175 0.863 0.766 0.784 0.715 1.03 0.915 
83 130 115 128 116 135 120 133 0.978 0.866 0.964 0.873 1.01 0.903 
66 278 228 282 234 283 230 235 1.18 0.971 1.20 0.996 1.20 0.980 
71 332 297 352 317 384 333 327 1.02 0.910 1.08 0.970 1.18 1.02 
42 56 41.5 60 44 65 46 58 0.965 0.716 1.03 0.759 1.12 0.793 
Gl 55 38 56 39 59 40.5 48 1.14 0.793 1.17 0.813 1.23 0.845 
BR-1 74 63 75 67.5 98 80 80 0.922 0.788 0.938 844 1.22 1.00 
BRA 304 267 331 208 390 40 890 0.780 0.686 0.850 0.765 1.00 0.872 
BR-5 57 50 68 61 92.5 4 76 80 0.713 0.625 0.850 0.763 1.16 0.950 
of the shell. The geometric ratios, yield points, experimental and shown in Fig. 2 for illustrative purposes. In constructing these 


theoretical collapse pressures for these cylinders are given in 
Table 2. The first six models in this table were tested some 20 
years ago by Windenburg and Trilling (11), although they did 
not investigate theoretically the effect of initial out-of-roundness 
on the collapse pressure. The last three models, which are multi- 
bay cylinders, have been tested recently at the Taylor Model 
Basin (13, 14). 

A comparison of columns 6 and 8 in Table 2 shows a consider- 
able discrepancy between the experimental and theoretical 
buckling pressures for even the simply supported cylinders. For 
some models, it would also appear that axisymmetric yielding 
rather than buckling was the controlling mode of failure. This 
can be seen by comparing columns 8 and 10in Table 2. However, 
it will be seen later when out-of-roundness is taken into account 
that the theoretical pressures for buckling-type failures are lower 
than the axisymmetric yield pressures. It is also of interest to 
note that Models 42 and 61 are the only models for which the ex- 
perimental buckling pressures are higher than those predicted 
theoretically for simply supported cylinders, although the same 
supports were used for these two models as for the four models 
preceding them. 

Using the geometric ratios and yield points shown in Table 2 
apd Equations [11] and [12] in Table 1, it is possible to construct 
curves showing the relation between the pressure at which yielding 
first occurs in the cylinder wall p,, and the ratio of the initial out- 
of-roundness to the shell thickness ¢/h. Two such curves are 


curves, the value of m used in Equation [12] was the value of 
m which minimized Equation {11]in Table 1. These values of m 
are listed in parentheses in columns 8 and 9 of Table 2, These 
values of m do not actually give the minimum p, for a given e/h. 
This point will be discussed later. It also should be noted that 
for e/h = 0 some of the curves in Fig. 2 do not attain the elastic 
buckling pressures tabulated in columns 8 and 9 of Table 2. 
When this is found to occur it means that the pressure to cause 
axisymmetric yielding of the shell is lower than the elastic 
buckling pressure. 

In Table 3 are tabulated the maximum e/h-values obtained 
using the different methods for determining the initial out-of- 
roundness described earlier. The e/h-values were determined at 
mid-length for the unstiffened cylinders and at mid-length of the 
bays for the multibay cylinders. Two values are listed for each 
model under methods (a) and (c) because, in these methods, the 
number of lobes into which the perfect cylinder would buckle is 
used and vais number is usually different for simply supported 
and clamped ends. 

Now, selecting the experimental buckling pressures listed in 
column 6 of Table 2 in conjunction with their corresponding ec- 
centricities listed in columns 2, 3, and 6 of Table 3, one can plot 
points on curves similar to Fig. 2 which represent values deter- 
mined experimentally. In Table 4 we also give a numerical com- 
parison of the theoretical and experimental pressures at which 
visible lobes first occur. The theoretical values in this table 
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TABLE 5 AVERAGE VALUES OF COLUMNS (6) Aa AND (8) IN 
TABLE 4, CLASSIFIED ACCORDING TO CYLINDER TYPE 


Cylinder ——Method (a)— —Method (b)— —Method (c)—~ 
type C.E. 8.38. C.E. 3.8. C.E. 3.3. 
Unstiffened....... 1.08 0.85 1.12 0.89 1.18 0.91 
Stiffened......... 0.85 0.75 0.88 0.79 1.08 0.93 


were obtained from the e/h-values listed in Table 3 in conjunction 
with the theoretical p, versus e/h curves similar to Fig. 2. The 
last three columns in Table 4 show the ratios of theoretically pre- 
dicted pressures for the occurrence of a visible lobe to those ob- 
tained experimentally, according to the. various methods used for 
determining the initial out-of-roundness. The average values of 
these last three columns, classified according to whether the 
cylinders were stiffened or not, are tabulated in Table 5. It can 
be seen from Tables 4 and 5 that use of method (a), with the 
assumption of simply supported ends, was the most conservative 
in most cases and predicted pressures which were always below 
those obtained experimentally. However, the best correlation 
between the experimental results and the simplified theories 
discussed in this paper appears to be obtained when method (c) is 
used for determining the initial out-of-roundness and the cylinders 
are assumed to be simply supported. 

It was mentioned earlier that for a given e/h the value of m that 
would give the lowest p, was not necessarily the value of m which 
minimized the expression for por (Equations [11] in Table 1). It 
was also noted at that time, however, that the error obtained by 
assuming this to be so was not very great and also greatly re- 
duced the computational labor. Some idea of the error involved 
can be obtained by referring to Fig. 3. The curves in this figure 
are plots of m which minimize e/hforagiven p,. Fora p, of 80 psi 
it can be seen that the minimum value of e/h is 0.34 at m = 18, 
while at m = 15 (the value which minimizes por) the value of 
e/h is 0.4. However, if we now fix the value of e/h at 0.4, then 
the minimum value of p, is 75 psi and occurs atm = 19. Similar 
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results are obtained with the other curves and also with the m 
versus e/h curves for clamped-end cylinders which we have not 
included here. It is thus seen that assuming the value of m which 
minimizes por will also minimize p, for a given e/h is slightly on the 
unsafe side, the magnitude of the error depending on the value of 
e/h taken. However, owing to the fact that actual cylinders under 
hydrostatic pressure collapse in substantially the same number of 
lobes as is predicted by theory for the perfect cylinder, and also for 
simplicity in calculations, we have ignored this discrepancy. 

It also will be remembered that in applying the analyses to 
stiffened cylinders we made the assumption that the stiffening 
rings were perfectly circular. This, of course, never occurs in 
practice. Some idea of the degree of circularity actually present 
can be obtained by reference to Fig. 4. This figure shows the 
initial circularity contours of the stiffening rings bordering, and the 
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shell at the center of one of the bays in which lobes first appeared 
in Model BR-4. If we adopt method (c) described earlier as our 
criterion for initial out-of-roundness, then the stiffening rings of 
this model had about one tenth the:out-of-roundness of the shell. 
Similar results also were obtained for the other models. Thus, for 
the problem of inter-ring collapse of stiffened cylinders, the as- 
sumption of zero initial out-of-roundness of the stiffening rings 
appears to be a reasonable one. 

It might be thought that a harmonic analysis of the initial out- 
of-roundness would show that the amplitude of the harmonic 
component corresponding to the value of m which minimized por 
(Equation [11] in Table 1) was many times that of the other com- 
ponents. To investigate this point a harmonic analysis, using 72 
subdivisions of the circumference, was made of the initial circu- 
larity contour at Station 4 of Model BR-4 using Filon’s method 
(15). The number of waves into which a perfect cylinder with 
shell dimensions similar to BR-4 would buckle is, according to the 
linear theory used herein, 11 for simply supported ends and 12 for 
clamped ends. From the harmonic analysis it was found that 
there were nine harmonics with amplitudes greater than the 
eleventh and two greater than the twelfth. It also was found that 
the amplitude of the largest harmonic (m = 8) was more than 
twice that of the twelfth mode and more than five times that of 
the eleventh mode. Further, even this largest amplitude was 
much too small to effect a reasonable correlation between theory 
and experiment. Thus the harmonic analysis of the BR-4 initial 
out-of-roundness contour did not produce any results which 
might be useful in practice. 

In conclusion, it might be of intereat to mention that some at- 
tempts have been made to determine experimentally the longi- 
tudinal form of the buckling displacemnts in the multibay cylin- 
ders. The results of these few investigations are summarized in 
reference (14). However, more experimental work still remains to 
be done before any conclusions regarding the shape of the buckling 
displacement can be made. 


SUMMARY 


In the preceding sections an attempt has been made to explain 
some of the discrepancies that exist between experimental and 
theoretical results for cylinders subjected to external hydrostatic 
pressure. To do this it was assumed that (a) the actual boundary 
conditions were somewhere between the extremes of simple sup- 
ports and clamped ends; (b) the initial out-of-roundness was 
similar in form to one of the modes into which a perfect cylinder 
would buckle; (c) the stress distribution in the equilibrium prob- 
lem for the perfect cylinder could be represented by the membrane 
stresses; and (d) the linear small-deformation equations of equi- 
librium would describe the problem adequately. Any cold-work- 
ing, residual, or welding stresses, or any elastic nonhomogeneity 
that might have been present were neglected. It also was as- 
sumed that failure (formation of a lobe) would occur when the 
stresses at the most highly stressed point satisfied the octahedral 
shear-stress criterion. Three simplified methods of measuring the 
initial out-of-roundness were also investigated. The simplified 
analyses, together with the different methods of measuring out-of- 
roundness, were applied to nine steel welded cylinders with length- 
diameter ratios of !/, to 2, thickness-diameter ratios of 0.0025 to 
0.0065, and yield points of the steel of 30,000 to 60,000 psi. The 
correlation between experimental and theoretical results was 
quite good, when method (c) was used for determining the initial 
out-of-roundness of the cylinders. 
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Appendix 


APPROXIMATE Equations OF EQUILIBRIUM FOR AN INITIALLY 
Out-oF-RounD CYLINDER 


These equations have been derived previously by Bodner and 
Berks (3) using a different approach. Our reason for this new 
derivation is to make clear what approximations and assumptions 
are involved in the simplified equations of equilibrium. 

Consider a cylinder of mean radius R which has a small radial 
initial out-of-roundness wo. Selecting z, 6, and z as co-ordinate 
axes and denoting the longitudinal, tangential, and radial elastic 
displacements by U, V, and W, it is then possible to show that 
the strains at the middle surface are given by 
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We shall simplify these complicated expressions to the following 
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For the changes in curvature, we shall use the quantities pre- 
viously used by Donnell (8) for perfect cylinders, viz. 
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The work done Wp by the uniform external pressure acting on 
all sides of the cylinder is given by the product of the pressure 
and the change in volume of the cylinder. Thus there follows 
(16) 
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We shall simplify this complicated expression for Wp to the follow- 
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It is also convenient to consider the U, V, and W displacements 
to be made up of two parts 
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where UO, V, and W are the displacements which would occur in 
the equilibrium problem of a perfectly circular cylinder under 
uniform external pressure. 

The total potential of the system U, is then obtained by add- 
ing the extensional and bending-strain energies of the shell and 
subtracting the work done by the external pressure. In calculat- 
ing the extensional energy, we retained the terms in U, Vv, W, Wo 
through the second order and the terms in U, 7, W (directly pro- 
portional to the applied pressure) through the first order only; 
also, we neglected the effect of the deflection of the shell between 
supports on the displacements and stresses of the perfect cylinder. 

Variation of Up with respect to u, v, and w then gives the dif- 
ferential equations of the problem. Further manipulation of these 
equations results in the following 
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where @,, dy, and-7,9 are the membrane stresses which would 
occur in a perfectly circular cylinder and which were assumed to 
be constant in deriving the foregoing equations. 

For the case of uniform external hydrostatic pressure applied 
on all sides of a perfectly circular cylinder, and neglecting the de- 
flection of the shell between supports, there results 
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Substituting the Relations [20] into Equation [19c] yields 
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Now define a stress function of the total membrane stresses 
F, such that 
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Using Equation [14] and the stress-strain relations it is then 
easy to show that 
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As we are interested here in the linear problem, and as we have 
previously neglected the deflection of the shell between supports, 
Equation [23] reduces to 


Equations [19a], [19b], [21), and [24] together with the appro- 
priate boundary conditions define the problem. 
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